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The momentum correlation between protons and lambda particles emitted from central Pb+Pb colli- 
sions at y^s^ =17.3 GeV was studied by the NA49 experiment at the CERN SPS. A clear enhance- 
ment is observed for small relative momenta (qi^y < 0.2 GeV). By fitting a theoretical model, which 
uses the strong interaction between the proton and the A in a given pair, to the measured data a value 
for the effective source size is deduced. Assuming a static Gaussian source distribution we derive an 
effective radius parameter of /?g = 3 -02 ± 0.20(stat.)l^;^(syst.) fm. 
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1 Introduction 

Heavy-ion collisions at ultra-relativistic energies produce strongly interacting matter under extreme con- 
ditions. The main goal is to create in these reactions a state in which the confinement of quarks and 
gluons inside hadrons is no longer effective, the so-called quark-gluon plasma. This strongly compressed 
matter undergoes a rapid expansion with a drop of temperature and energy density. Two-particle momen- 
tum correlations provide unique information on the size and dynamic evolution of this fireball and are 
therefore a widely employed observable in heavy-ion physics. Usually correlations of identical charged 
pions are studied which, due to the high available statistics, allow a multi-dimensional study of radius 
parameters |IT]|2]|3l. Less frequently two-proton correlations are analyzed E] |5] ID |7l, although only 
one-dimensionally. Moreover, the large abundance of strange particles produced in heavy-ion collisions 
allows to study also two particle correlations between pairs of strange particles or pairs of strange and 
non-strange particles. For example the correlations of identical kaons were investigated at the CERN- 
SPS [8] as well as at RHIC [9J. In this paper we report on the measurement of the p — A correlation 
function in central Pb-i-Pb colhsions at ^/s^ - 17.3 GeV at the CERN-SPS. 

It was suggested that also the momentum correlation between A and protons can be employed to measure 
the size of the emitting source lITOl . The correlation function of p — A pairs is only affected by the strong 
interaction between the particles. This distinguishes p — A correlations from the two-proton case, for 
which the correlation function is dominated by the repulsive Coulomb interaction and the Fermi-Dirac 
statistics at low relative momenta. Both effects are absent in the p — A correlation function, which should 
therefore be more sensitive to large source sizes [|TOl . However, the knowledge of the strong interaction 
between protons and A is necessary to relate the strength of the correlation to the size of the emitting 
source. There is a substantial set of data available on low-energy elastic A p scattering nTl [T2l[T3l . and 
on K^d — )■ Ap %^ [ 141 [15] |T6J, as well as pp —5- pK+A 1 17, 1 81 reactions. Also A hypemuclei provide 
important information on the A nucleon interaction. Based on this data many theoretical analyses derived 
p — A scattering lengths and effective interaction ranges ||T9]l20ll2ni22ll23ll24ll25l which can be used to 
calculate the p — A correlation function. 

A preliminary study by the NA49 experiment at the SPS was reported in f26l. Here we describe the final 
results of a new and improved analysis L27 J . Similar studies of p — A correlations in heavy ion collisions 
were performed at lower [|28l |7] and higher ||29l center-of-mass energies, allowing to investigate the 
evolution of p — A correlations with JJ^. 



2 Data Analysis 

The data presented here were measured by the NA49 experiment at the CERN SPS. A detailed description 
of the experimental setup can be found in ll30l . Charged particles produced by interactions of the Pb 
beam in a thin Pb-foil target are tracked with four large- volume Time Projection Chambers (TPCs). Two 
TPCs are placed inside two superconducting dipole magnets, while the other two are situated outside 
of the magnetic field. Since the latter measure long pieces of the particle tracks, they allow a precise 
determination of the specific energy loss dE /Ax inside the detector gas (typical resolution of 4 %) and 
thus particle identification in a large region of phase-space. Additional particle identification is provided 
around mid-rapidity by Time-Of-Flight detectors. A Zero Degree Calorimeter is used to measure the 
energy in the projectile fragmentation region from which the centrality of the reaction can be deduced. 
This analysis is based on 2.8 • 10^ Pb-i-Pb events at ^/s^ = 17.3 GeV recorded in the year 2000, which 
cover the 23.5 % most central part of the total inelastic cross section, corresponding to an averaged 
number of wounded nucleons of (A'w) = 262. 
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Fig. 1: The invariant mass distribution of A candidates in two exemplary phase-space bins ((a) — 1.2 < y < —0.4 
and 0.6 <pt<1.2 GeV/c, (b) -0.4 < y < 0.4 and 1.8 < pt< 2.4 GeV/c) for central Pb+Pb reactions at = 
17.3 GeV. The two vertical lines indicate the mass windows used to define A candidates. 



2.1 A reconstruction 



The A hyperons are detected via their charged decay A — )• pn^, using the same methods as described in 
Il3ni32l . The A reconstruction is done by forming pairs of positively and negatively charged tracks and 
extrapolating them towards the main interaction vertex. The positively (negatively) charged tracks are re- 
quired to have at least 50 (30) reconstructed points. Pairs with a distance of closest approach of less than 
0.5 cm anywhere between the position of the first measured point on the tracks and the target plane are 
considered as V° candidates. Assigning proton and pion masses to the positively and negatively charged 
decay particle, the invariant mass of a A candidate is calculated. A significant reduction of the combi- 
natorial background can be achieved by applying several selection criteria to the A candidates. In this 
analysis it is required that the secondary vertex is separated by at least 25 cm in beam-(z)-direction from 
the target plane. Additionally, the back-extrapolation of the flight path of the A candidate must not devi- 
ate from the interaction vertex position in the transverse directions x and y by more than \Ax\ = 0.75 cm 
and \Ay\ = 0.375 cm. The signal-to-background ratio is further improved by enriching the protons in 
the sample of positively charged tracks by applying a momentum-dependent cut on the measured energy 
loss (dE/dx). An additional dE/dx cut on the negatively-charged tracks also allows to reject electrons 
from photon conversions. 

An important point with respect to correlation studies is the requirement that each A candidate must be 
unique. If it happens that a daughter track of a given A candidate is also assigned to another one, a 
strong artificial correlation between both candidates is created which in turn affects the measured p — A 
correlation function. In this analysis it is therefore ensured that any given track is used only once as a 
daughter track. Similarly, A daughters that were not reconstructed as a single track, but as two track 
pieces (split tracks), will cause a distortion of the measured correlation function. To exclude these tracks 
it is required that the number of measured points of each accepted track is higher than 50 % of the number 
of points that this track could maximally have according to its trajectory in the TPCs. 

Figure [Ushows the distribution of invariant mass m,„,, of A candidates obtained after assigning the proton 
(71^) mass to the positive (negative) daughter track for two intervals of the center-of-mass rapidity y 
and transverse momentum pt. A are accepted in a mass window [niQ — Am, mo + Am] of a half-width of 
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Fig. 2: The purity PA{y,Pt) of the selected A candidates as function of rapidity y and transverse momentum pt. 

Am = 2 MeV/c^, where niQ = 1.115683 GeV/c^ is the literature value for the A mass Il33ll . This mass 
window was chosen in order to optimize the signal-to-background ratio and reduce the corrections for 
the signal purity. 

The ratio between signal and background inside this mass window and with it the purity of the A sample 
vary over phase-space. A ;t^-fit to the invariant mass distributions in different y and pi bins was per- 
formed in order to determine the relative contributions. The fit function used is the sum of a polynomial 
for the combinatorial background and a function for the A signal. The measured shape of the A signal 
results from a convolution of the resolutions for the single track momenta and the secondary vertex posi- 
tions and is found to be well described by an asymmetric Lorentz-curve. By subtracting the background 
function from the measured invariant mass distribution the number of real A is determined. Defining 
the purity Pa by the ratio of the number of real A to the number of accepted A candidates in the chosen 
mass window one obtains the result shown in Fig. |2] The phase-space averaged value of the A purity is 
(Pa) =69%. 



2.2 Proton reconstruction 

Protons are identified via their energy loss dE/dx as measured in the two Main-TPCs. A valid proton 
track is required to have at least 50 reconstructed points. An additional cut on the impact parameters 
in the target plane {\b-^\ < 5.0 cm, \by\ < 2.0 cm) reduces the contribution from secondary tracks. In 
order to assign the probability of being a proton to a given track, the energy loss spectra measured in 
bins of total momentum p are fitted by a sum E{x,p) of asymmetric Gaussians using a minimizing 
procedure ll34]| : 



1 



ni 



exp 
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2 V(l±5)a,-,,(p) 



(1) 



Here, Ai{p) denotes the yield for particle type /, ni the number of tracks in a given track length interval 
/, Xi{p) the most probable dE/dx values for particle type /, x the measured dE/dx value of the track 
under consideration, Gij{p) the width of the Gaussian, and 5 the asymmetry parameter. The parameters 
Aj{p), Xj{p), and (yj[{p) are determined by the fitting procedure in each p bin separately. The widths for 
the other particles types and different track length bins are derived from the parameter Onip) for pions: 
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Fig. 3: The AE /&>c spectra in two exemplary total momentum bins ((a) 4.0 < p < 5.0 GeV/c, (b) 12.6 < p < 
15.9 GeV/c) together with their decomposition into contributions from protons, kaons, pions, and electrons (from 
left to right). Shown is the measured energy loss normaUzed to the minimum ionizing value dE/diMip- The 
vertical lines represent the upper dE/dx cuts that result in a purity of Pp = 0.8. 

'^i.iip) = (^n{p) ■ (■^(/■^7t)"(l/ Vl)- The exponent was determined to be a = 0.625 and the parameter d is 
fixed to 0.065. Results of the fit procedure for two different total momentum bins are shown in Fig. [3l 
Since for low particle momenta the energy loss curves of different particle species cross each other and 
particle identification is thus not possible, a lower cut on the total momentum of p > 4 GeV/c is applied. 
To exclude the region of the Fermi plateau also momenta above 50 GeV/c are discarded. Based on 
these fits a momentum dependent cut on the measured dE/dx values for single tracks is defined such that 
the accepted tracks always have the same probability of being a proton. For the standard analysis this 
probability is set to 80 %, equivalent to a constant proton purity Pp = 0.8. As in the case of the A decay 
daughters it is ensured that split tracks are removed from the track sample by rejecting tracks that have 
less than 55 % of the number of geometrically possible points. 

2.3 Determination of the p — A correlation function 

The selected A and proton candidates are then combined to form p — A pairs. In order to avoid trivial 
auto-correlations a track that is used to reconstruct the A candidate is removed from the primary proton 
sample. The pair distribution ^(^inv) is measured as a function of the generalized invariant relative 
momentum of the p— A pair ^inv, which is defined as the modulus of ^ = q — P{qP)/P^, with q = pp — p^, 
P = Pp+ Pa, and qP = — m\. Here, pp and p^ are the 4-momenta of the proton and the A. In the two 
particle center-of-mass system q reduces to {0,2^*}, with 2k* being the 3-momentum difference in this 
reference frame ll35l . 

An event-mixing method that combines proton and A candidates taken from different events, is employed 
for the construction of the uncorrelated background B(ginv)- The measured correlation function is thus 
defined as: 



The normalization constant A'^ is determined by requiring Cmeas(^inv) = 1 in the region 0.2 < < 
0.3 GeV. Since the reconstruction of real pairs is affected by the limited two-track resolution of the 
detector 1^3, a distance cut between the track of the primary proton and the track of the positive A decay 
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Fig. 4: The phase-space population of protons (a), A (b) as a function of rapidity y and transverse momentum 
pt. The distribution for the resulting p — A pairs (c) is shown versus pair rapidity = j In gptgA-pz'p-Pz A 



particle is applied for real pairs as well as for the mixed event pairs. For each pair it is required that the 
tracks have an average separation of at least 3.0 cm. This average is determined as the arithmetic mean 
of the track distances determined in planes perpendicular to the beam axis. For each TPC two planes 
are taken into account. Their distances to the target plane are 160.5 cm, 240.5 cm, 540.0 cm, 620.0 cm, 
910.0 cm, and 990.0 cm. Pairs that do not pass the distance cut are discarded. 

Figure |4] shows the phase-space population of the accepted protons. A, and p — A pairs. Averaging 
over all measured p — Apairs we find {kt) = 0.53 GeV/c with = ^|a,p + A,a| {{m^ = 1.18 GeV with 

mt = Y^^f + (i(mA + mp))2). 

The measured correlation function can be affected by the finite momentum resolution of the detector. 
In 131 an extensive investigation of its influence on the radius parameters extracted from correlations 
of identical charged pions as measured with the NA49 experiment is discussed. Due to the excellent 
momentum resolution of the NA49 detector, it turned out that the impact is negligible and a correction 
for this effect is not necessary. Even though the momentum resolution for a A ((cJ^^) < 1%) is worse 
than for the primary track, the resulting effect on the measured p — A correlation function is still clearly 
smaller than all other systematics effects. Therefore no correction is applied in this analysis. 

A substantial fraction of the measured protons and A originate from weak and electro-magnetic decays 
of heavier particles (feed-down). In the following it will be assumed that these decay particles are 
not correlated, because the decays happen long after the thermal freeze-out, and will thus reduce the 
observed p — A correlation function. For the A the feed-down originates from S^, S*', and Z" decays, 
while in case of the protons the decays of A, and £+ contribute. We calculate the fraction of protons 
and A originating from feed-down {F^{y,pi) and FaCJjA)) via a simulation procedure. A, E^, and 
are generated according to their measured phase-space distributions fST]. The measured A include the 
which cannot be separated experimentally. For the E" the same input distributions are assumed as 
for the E^ scaled by the ratio E^/E^ of the total multiplicities taken from statistical model fits [36|. 
The daughter tracks of the generated particles are followed through the NA49 detector setup using the 
Geant3.21 package |[37ll . The response of the TPCs to the traversing particles is simulated with NA49 
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Fig. 5: The fraction of protons iFp{y,pt)) (a) and of A (FA{y,pt)) (b) originating from feed-down as a function of 
y and pt. 



specific software. In a next step the simulated raw signals are added to measured raw data and processed 
by the same reconstruction program as used for the experimental data. By applying the same cuts as in 
the normal analysis, the feed-down contributions to the measured A and protons are determined. The 
y and pt dependences of the feed-down fractions Fp{y,pi) and Fx{y,Pt) are summarized in Fig. |5] The 
phase-space averaged values are (Fp) = 22 % and {F^) = 43 %. 

The final p — A correlation function Ccon-(9inv) results from the measured Cnieas('?inv) by applying a com- 
bined correction factor {K{q\av)) for purity and feed-down. This factor is determined by averaging the 

product Pp(3;p,;?t,p)^A(3'A,A,A) (1 -^p(3'p, A,p)) (1 -^a(3'a,A,a)) over all reconstructed p - A pair com- 
binations falling into a given bin of relative momentum ^inv- The corrected correlation function thus 
follows from: 

C'meas('?mv) 1 



Ccorr(9in 



(^(^inv)> 



+ 1. 



(3) 



3 Results 

In total 70920 p — A pair candidates, corresponding to 17520 real p — A pairs (i.e. after correcting 
for purities and feed-down), with ^inv < 0.2 GeV were measured. Dividing the signal distribution by 
the event-mixing background and correcting with the purities and feed-down contribution according to 
Eqs. ^ and Q yields the final p — A correlation function as shown in Fig. |6] The correlation function 
exhibits a significant enhancement for small ^inv- Such a correlation would be expected as the effect of 
the strong interaction between the proton and the A. 

3.1 Fit to theoretical calculation 

Since the shape of the momentum correlation function was shown by Wang and Pratt to depend on 
the size of the emitting source lITOl . it can be used to extract its radius. The necessary prerequisite is 
a quantitative knowledge of the p — A interaction. Here we use a functional form of the theoretical 
correlation function Qh that is based on the model of Lednicky and Lyuboshitz |[38l[35l . It employs an 
effective range approximation of the S-wave p — A interaction. The source size is required to be larger 
than the effective range of the interaction. The strength of the interaction is defined by four parameters: 
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Fig. 6: The corrected p — A correlation function for central Pb+Pb reactions at y/s^ = 17.3 GeV, shown as a 
function of the invariant relative momentum q^^y. The data represent an average over the whole acceptance of the 
NA49 experiment. The lines display results of the fit with a theoretical correlation function 1^ 38] |3?1 (see text for 
details). Only statistical eiTors are shown. 

the effective ranges (cIq) and the scattering lengths /q (/,}) for the singlet S = (triplet 5=1) state. 
In our fits we use values of = 2.92 fm, d^ = 3.78 fm, /q = -2.88 fm, and = -1.8 fm, as suggested 
in lITOll . Under the assumption of unpolarized particle production, the relative contribution of pairs in the 
singlet and the triplet state is 1 : 3. Furthermore, we use the same spherically symmetric Gaussian spatial 
distribution S{r), for both the proton and the A source: 



Then, the theoretical correlation function Qh assumes the functional form as quoted in reference f29]. 
By fitting it to the data one obtains the effective radius parameter Rq and an additional parameter A that 
takes possible reductions of the height of the correlation into account: 



These reductions occur, if the correction for the particle purities and the feed-down, as defined in 
Eq. (|3]l, are insufficient. If both parameters are left to vary freely, the best fit is obtained for Rq = 
2.70ib0.60(stat.) fm and A = 0.77 ib0.38(stat.) (green dashed line in Fig.[6ll. The agreement within er- 
rors of the fitted A value with unity underlines the consistency of the correction procedure. This justifies 
to fix A = 1 and to use only Rq as a free parameter, which reduces the resulting error. With this constraint 
we obtain Rq = 3.02 ± 0.20(stat.) fm (blue sohd hne in Fig. O. 

3.2 Systematic uncertainties 

Systematic errors arise both in the extraction of the correlation function and in the theoretical model. 
The former uncertainties were studied by making small changes in the analysis procedure. By varying 
the corresponding cuts on the measured energy loss used to identify the primary protons, as well as 




(4) 



Cfit(<?inv) = A (Cth(<?inv) - 1) + 1. 



(5) 
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Fig. 7: The corrected p — A correlation function CcorrC^inv) for central Pb+Pb reactions at y/s^= 17.3 GeV shown 
as a function of the invariant relative momentum qi^y for different purities of protons (a) and A (b), as well as 
different feed-down contributions to protons (c) and A (d) (see text). The lines display the results of fits with a 
theoretical correlation function ||38] l35l . 

the decay protons of the A candidates, the particle purities Pp{pp) and PAiyAiPtA) can be varied to a 
certain extent. This changes the measured correlation function Cmeasi^mv)- However, after applying 
the appropriate correction factor (7c'(^inv)) the same corrected correlation function should be obtained. 
Figure |7] shows comparisons of Ccorr(^inv) for different proton (Fig. |7]a)) and A purities (Fig. |7jb)). 
Even though the correction factor changes quite dramatically (?a 45 % in case of the proton purity), the 
resulting correlation functions agree quite well. The systematic error on the effective radius parameter 
Rq is derived by taking the maximal difference of Rq as obtained by fits to the different correlation 
functions. It is found that it changes between 2.91 — 3.28 fm (2.98 — 3.29 fm), if the proton (A) purity is 
varied between 50 - 90 % (64 - 72 %). 

A similar study is performed by varying the feed-down contribution by changing the cuts on the impact 
parameter bx and by (Ax and Ay) of the proton tracks (Fig. |7]c) and A candidates (Fig. |7jd). One finds 
that Rq changes maximally between 3.02 — 3.13 fm (2.92 — 3.23 fm), if the contributions to the protons 
(A) from feed-down is varied between 17 — 25 % (39 — 46 %). 

The dependence of the radius parameter on the region in ^inv that is used to determine the normalization 
constant (see Eq. ([2])) is investigated by varying the size and position of this region. It is found that Rq 
changes by maximally ±0. 1 1 fm. 

By taking the quadratic sum of the different contributions a total systematic error on Rq of +0.44 fm and 
— 0.16 fm is estimated. 
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Fig. 8: The effective radius parameter Rq for p — A correlations as a function of center of mass energy y/s^. The 
figure also includes data on Ar+KCl collisions by the HADES collaboration [28 1 and on central Au+Au collisions 
by the E895 Q and the STAR 1291 experiments. The systematic errors are represented by the brackets. Please 
note that the STAR result corresponds to a slightly higher pair (mt) than for the NA49 measurement. Assuming 
the nil dependence shown Fig. |9]the STAR data point would move up by 6 %, if it was measured at the same 
(mt). 

Another source of systematic uncertainty of the extracted radius parameter arises from the limited pre- 
cision of the knowledge of the scattering lengths and effective ranges used in the calculation of the 
theoretical correlation function. Therefore, the fits were repeated with dg, jj, /q , and taken from 
|[ni[T8l|20l|25l|22l|23l[33. The largest deviation from the Rq value obtained with the standard parame- 
ter set [10] is observed when using the parameters extracted by the COSY-11 collaboration lITTl . While 
in all other cases the difference is smaller than ±0.1 fm, it is found to be +0.274 fm for the COSY- 
1 1 parameters which is still close to our statistical error. Therefore we conclude that the choice of the 
parameters describing the p — A interaction has a negligible effect on the final result. 

3.3 Energy dependence of the effective radius parameter 

The result of this study is compared to data at lower and higher center-of-mass energies in Fig. [8] Good 
agreement with the effective radius parameters measured for central Au-i-Au collisions by the E895 col- 
laboration at y/s^ - 3.83 GeV [7] and by the STAR experiment at ^/s^ = 200 GeV Ii29,l is observed, 
indicating that there is no significant change from AGS to RHIC energies. A similar observation was 
made for correlations of identical charged pions |j3]- 

At even lower energies (y^^ = 2.61 GeV) the HADES collaboration measured a significantly smaller 
effective radius parameter for p — A correlations of Rq = 2.09 ±0. 16(stat.) fmf28l. However, the fireball 
volume in the case of Ar-i-KCl collisions is expected to be smaller than for Pb-i-Pb and Au-i-Au collisions. 
In fact, it was demonstrated in |[28l that the measured p — A radius parameter is dominated by the re- 
action geometry and scales approximately as A^/^, irrespective of center-of-mass energy. A comparable 
observation was previously made for two-proton correlations in the target fragmentation region ||6]. 
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Fig. 9: The effective radius parameter Rq extracted from the correlation functions of n^n^ 13], and K^K^ 

10. PP Hi and p — A pairs versus (mt) as measured by NA49 for central Pb+Pb collisions at ^/s^ = 17.3 GeV. 
The systematic errors are represented by the brackets. For explanation of curves see text. 

3.4 (mt) dependence of the effective radius parameter 

Figure |9]shows a comparison of Rq, as determined in this analysis from the p — A correlation function, to 
radius parameters derived from correlations of tt^ pairs [3| and charged kaon pairs |'8l, as well as from 
two-proton correlations [4 |, at different average transverse masses {m^). In case of the n^n^ and KK 
correlations, Rq was calculated from the three-dimensional radius parameter components as: 



The decrease of the effective radius parameters with increasing transverse mass is generally attributed 
to the presence of collective flow in the fireball. The measurement of the p — A correlation allows to 
extend this study to higher (mt) . In fact, Rq for p — A pairs is significantly smaller than the effective 
radius parameter extracted for pions and kaons at lower (mi) and is thus in agreement with the expected 
behavior. The dashed curve in Fig. |9] corresponds to a simple oc (wit)^'/^ dependence. The solid line is 
based on the following mt dependences of the three radius components, as suggested by hydrodynamical 
approaches ||39l |40] 1411 : 



Here Rg^o is the transverse size of the particle source, j8t = Vf/c the transverse flow velocity, T]f the 
transverse flow rapidity Tjf = (1/2) log(l + Vf)/(1 — Vf), T the kinetic freeze-out temperature, Tq the total 
lifetime of the source and At the emission duration. Under the assumption that all particle species freeze 
out from the same expanding source, we calculate Rq from Eq. Q, using Eqs. (|7]l-(|9ll. The solid curve in 
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Fig. |9] corresponds to the parameter T]f = 0.8, Tq = 0.6 fm/c, At = 3.4 fm/c and T = 90 MeV, as extracted 
by fits with a blast wave model to the pion correlations [3J. The overall normalization has been adjusted 
to fit the data. A reasonable description of the effective radius parameters for most particle species can 
thus be achieved, with the notable exception of the two-proton correlation. 

4 Summary 

We report on the measurement of the p — A correlation function in momentum space for central Pb+Pb 
collisions at y^^J^ = 17.3 GeV. The p — A pairs exhibit a clear positive correlation for small relative 
momenta. By comparison to a calculated correlation function a one-dimensional Gaussian source size 
parameter of Rq = 3.02ib0.20(stat.)lQ |g(syst.) fm is determined. This value is in good agreement with 
measurements for Au-i-Au collisions at lower and higher center-of-mass energies. The (mt) dependence 
of the effective p — A radius parameter follows the expectation for an expanding source as described by 
hydrodynamics, when compared to other two-particle correlation results. 
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